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1. Introduction 

The discovery and identification of new degrees of freedom and interactions at high-energy 
colliders relies on the detailed understanding of Standard Model (SM) background pro- 
cesses. Prominent among these is the production of electroweak bosons (W, Z) in associa- 
tion with jets, of which one or more possibly contain bottom quarks. The prime example 
is the observation of the top quark at the Fermilab Tevatron collider, produced either in 
pairs Ijl], H or singly ^, since both of these mechanisms typically lead to W plus 6-jets 
signatures. Other crucial examples involve the search for a SM Higgs in association with 
vector bosons iWH/ZH), with the subsequent Higgs decay into a hb pair, a sought-for 
discovery channel both at the Tevatron Q and at the LHC [P, . Finally, in models which 
feature an extended Higgs sector, such as the MSSM or more generally a two-Higgs dou- 
blet model, a typical Higgs discovery channel is through Hhb and Abb final states, with 
an H/A — J- r^r^ decay. In this case, the SM process £~^i~bb can provide an important 
reference measurement. 

Next-to-leading-order (NLO) QCD calculations for the production of a vector boson 
in association with jets have by now quite a successful record. Accurate predictions for W 
plus up to four light jets 11,1^, 13| and for Z plus up to three light jets 15| 



have become available in the past few years. Associated production with heavy quarks, 
and in particular with bottom quarks, has been studied using various approximations. The 
Wbb and Zbb processes have been calculated at the NLO for the first time in refs. |16| 



and 1 17 1 respectively, by setting the bottom-quark mass equal to zero. Such calculations 
can be used only for observables that contain at least two 6-jets. The same processes 
have been considered again in refs. pi, E0|, where a non-zero bottom-quark mass has 



been used; however, the matrix elements still involved on-shell vector bosons, thereby 
neglecting spin correlations of the leptons emerging from W and Z decays. In the case of 



W production, this limitation has been recently lifted in ref. |21], which presents the NLO 
calculation for the leptonic process iubb. Other NLO calculations for final states with one 
6-jet, Wb and Zb and one 6-jet plus a light jet, Wbj and Zbj |2^, are also 
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available in the five- flavour scheme. All such calculations have played a role and/or have 
been extensively compared to the data collected at the Tevatron |26, 28, and now 



start to be considered in LHC analyses as well [30, 31]. 

In this paper we present a calculation of ivbb production that includes NLO QCD 
corrections (analogous to that of ref. [21|), and the first calculation at the NLO of i^i^bb 
production with massive bottom quarks; we retain the full spin correlations of the final- 
state leptons^. Furthermore, we match both of these results to the Herwig event generator 
by adopting the MC@NLO formalism |3^. Therefore, our results include all the relevant 
features which are important in experimental analyses, and can be used in order to obtain 
NLO predictions for a large class of observables, including those with zero, one and two 
6-jets. All aspects of the calculations are fully automated and analogous to the calculation 
recently appeared for Hti/Ati production ||3^. One-loop amplitudes are evaluated with 
MadLoop []3^ , whose core is the OPP integrand reduction method |35| as implemented 
in CutTools |3(:]. Real contributions and the corresponding phase-space subtractions, 
achieved by means of the FKS formalism |37|, as well as their combination with the one- 
loop and Born results and their subsequent integration, are performed by MadFKS p8| . 
The MC@NLO matching is also fully automated, and allows us to simulate for the first 
time ivbb and i^i~bb production with NLO accuracy, including exactly spin correlations, 
off-shell and interference effects, and hadron-level final states. All the computations are 
integrated in a single software framework, which we have dubbed aMC@NLO in ref. p3|| . 
We point out that the on-shell-Ty result of ref. [^] has recently been matched to showers 
in ref. || in the framework of the POWHEG box |0|. 

The phenomenology of ii/bb and i^i^bb final states is very rich and in fact transcends 
the prosaic role of background for Higgs or top-quark physics. Thanks to the aMC@NLO 
implementation, several QCD issues interesting on their own can now be addressed theoret- 
ically and the results efficiently compared to experiments. In this work we limit ourselves 
to providing new evidence that reliable and flexible predictions for the ii'bb and i^i~bb 
processes need to include: 



• NLO corrections; 

• bottom quark mass effects; 

• spin-correlation and off-shell effects; 

• showering and hadronisation. 



Detailed studies of these processes as backgrounds to specific signals, such as single-top 
and Hbb or Abb production respectively, are left to forthcoming investigations. 

This paper is organised as follows. In the next section we present several distributions 
relevant to iubb and i^i~bb production at the LHC, and report the results for total rates 
at both the Tevatron and the LHC. By working with a non-zero bottom mass, we are able 

^In the rest of this paper, ivbb and t^l^hh production as predicted by our simulations may also be 
denoted by Wbh and Zbb respectively. 
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to obtain predictions for the cases in which one or two 6's are not observed, and can thus 
have arbitrarily smah transverse momenta. We also give one example of the comparisons, 
at the level of hadronic final states, between the HW and HZ signals and their respective 
irreducible backgrounds which we have computed in this paper. We draw our conclusions 
in sect. ^. 

2. Results 

At the leading order (LO) in QCD lubb and i^i^bb production at hadron colliders proceed 
through different channels. Both final states can be obtained via a Drell-Yan-type mecha- 
nism, i.e., qq^'^ annihilation in association with a gluon splitting in a bb pair, see fig. |l|(a). 
Zbb, however, can also be produced by gluon fusion, see fig. ||(b), a channel that at the 
LO contributes a 30% of the total rate at the Tevatron, but turns out to be the dominant 
one (80%) at the LHC, owing to the larger gluon luminosity there. As we shall see in the 
following, the fact that Wbb and Zbb production are dominated by different channels at 
the LHC leads to important differences in the kinematical properties of final states, and in 
particular of 6-jets. 

We start by presenting results for the total cross sections at both the Tevatron, ^/s = 
1.96 TeV, and the LHC, = 7 TeV; the ivbb results are the sums of the i'^vbb and l^vbb 
ones (due to virtual- M/^"'' and W~ production respectively). In our computations we have 
set the lepton masses equal to zero, and is therefore not necessary to specify their flavour, 
which we generically denote by £ (for the charged leptons) and v (for the neutrinos); we 
always quote results for one flavour. For the numerical analysis we have chosen: 

with W = iv and U' = £^i~ in the case of Wbb and Zbb production respectively; the 
value of the 6-quark mass is that of the pole mass, = 4.5 GeV. We have used LO and 
NLO MSTW2008 four-flavour parton distribution functions for the corresponding cross 
sections, and the SM-parameter settings can be found in table ||. Given that the primary 
aim of this paper is not that of presenting precise comparisons with data, but rather 



Parameter 


value 


Parameter 


value 


mz 


91.118 




132.50698 


mw 


80.419 


Gf 


1.16639-10-^ 


rub 


4.5 


CKMij 




rat 


172.5 


Tz 


2.4414 




0.133551 


Tw 


2.0476 


(NLO,4) f N 


0.114904 







Table 1: Settings of physical parameters used in this work, with dimensionful quantities given in 
GeV. 
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Figure 1: Representative diagrams contributing to ii/bb and i^£~bb production at the leading 
order, (vbb production can proceed only via a qq' channel, diagram (a). For i'^i^bb production the 
qq channel, diagram (a), is dominant at the Tevatron, while the gg channel, diagram (b), largely 
dominates at the LHC. 





Cross section (pb) 


Tevatron ^=1.96 TeV 
LO NLO K factor 


LHC =7 TeV 
LO NLO K factor 


iubb 
£+e-bb 


4.63 8.04 1.74 
0.860 1.509 1.75 


19.4 38.9 2.01 
9.66 16.1 1.67 



Table 2: Total cross sections for ivbh and £~^£~bb production at the Tevatron {pp collisions at 
y/s = 1.96 TeV) and the LHC {pp collisions at ^/s = 7 TeV), to LO and NLO accuracy. These 
rates are relevant to one lepton flavour, and the results for liybb production are the sums of those 
for t^vbb and (r Dbb production. The integration uncertainty is always well below 1%. 

that of studying the defining features of the production mechanisms, in the CKM matrix 
(relevant to the Wbb results) we have neglected off-diagonal terms: this cannot change the 
conclusions we shall arrive at, but helps reduce the computing time. It should be clear that 
this is not a limitation of the code, since a non-diagonal CKM matrix can simply be given in 
input if one so wishes. Our runs are fully inclusive and no cuts are applied at the generation 
level, except for rrig+f- > 30 GeV in the i^i~bb sample. The predicted production rates at 
the Tevatron and at the LHC are given in table ^ where, for ease of reading, we also show 
the fully inclusive K factors. The contribution of the gg — )• Zbb + X channels is clearly 
visible in these results: at the Tevatron a{£~^i~bb)/a{£i'bb) is quite small (and of the same 
order of the ratio of the fully-inclusive cross sections a{Z)/a{W)), whereas at the LHC 
£^£~bb and £iybb differ only by a factor of two. 

We now study the impact of NLO QCD corrections on differential distributions, at 
both the parton level and after showering and hadronisation, and in doing so we limit 
ourselves to the case of the LHC, where the kinematical differences between Wbb and Zbb 
production are more evident. The parton shower in aMC@NLO has been performed with 
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Figure 2: Fractions of events (in percent) tliat contain: zero 5-jets, exactly one 6-jet, and exactly 
two 6-jets. The rightmost bin displays the fraction of 6-jets which are &&-jets. The two insets show 
the ratio of the aMC@NLO results over the corresponding NLO (solid), aMC@LO (dashed), and 
LO (symbols) ones, separately for Wbb (upper inset) and Zbb (lower inset) production. 



fortran Herwig [||, version 6.520^. 

We start by summarizing our results for 5-jet rates. Jets are reconstructed at the parti- 
cle level. In the case of MC simulations, this means giving all final-state stable hadrons^ in 
input to the jet algorithm. We adopt the anti-fer jet clustering algorithm [49| with R = 0.5, 
and require each jet to have PrU) > 20 GeV and \r]{j)\ < 2.5. A 5-jet is then defined as a 
jet that contains at least one 6-hadron; a 66-jet is a jet that contains at least two 6-hadrons 
(hence, a 66-jet is also a 6-jet). This implies that we make no distinction between the b 
quark and antiquark contents of a jet. We point out that at least another definition of 
6-jets exists |^0| which has a better behaviour in the — )• limit, in the sense that it 
gives (IR-safe) results consistent with the naive picture of "quark" and "gluon" jets. In 
practice, this is relevant only in the Pt 3> rrib limit. Since this region is not our primary 
interest in this paper, we stick to the usual definition; however, it should be obvious that 
any jet definition can be used in our framework. 

In fig. m we present 6-jet rates, as the fractions of events that contain zero, exactly 
one, or exactly two 6-jet(s). In the case of MC-based simulations, there are also events 
with more than two 6-jets and more than one 66-jet, but they give a relative contribution 
to the total rate equal to about 0.4% (for Wbb) and 0.6% (for Zbb), and are therefore not 
reported here. The rightmost bin of fig. |2| shows the fraction of 6-jets which are 66-jets. 
There is an inset for each of the two histograms shown in the upper part of fig. Each 
of the insets presents three curves, obtained by computing the ratio of the aMC@NLO 



^Automation of the matching to parton shower in the MCQNLO formalism to Herwig++ and to 
Pythia (see refs. and respectively) is currently under way. 

^In order to simplify the Herwig analyses, weakly-decaying B hadrons are set stable. 
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results over the NLO (solid), aMC@LO^ (dashed), and LO (symbols) corresponding ones. 
The 6-jet fractions are fairly similar for Wbb and Zbb production, and the effects of the 
NLO corrections are consistent with the fully-inclusive K factors. On the other hand, the 
66-jet contribution to the 6-jet rate is seen to be more than three times larger for ii'bb 
than for i^i^bb final states. This fact is again related to the different mechanisms for the 
production of a bb pair in the two processes considered here. At variance with the case of 
iiybb production, in a i~^l~bb final state the two 6's may come from the separate branchings 
of two initial-state gluons, and thus the probability of them ending in the same jet is much 
smaller than in the case of a g ^ bb final-state branching. We conclude this discussion 
by pointing out that the zero- and one-6-jet rates can only be obtained with a non-zero 
6-quark mass, since the one or two "untagged" 6's must be integrated down to pr = 0, and 
hence mf, 7^ is required in order to screen initial-state collinear divergences. This fact 
is a severe test condition for the computer programs used in the computations, because it 
may induce numerical instabilities. We stress that we did not impose low-px cuts on any 
of the final-state particles in MadLoop, MadFKS, and in the generation of hard events 
in aMC@NLO; our results are therefore completely unbiased, which is what gives us the 
possibility of computing quantities such as those reported in table ^ and in fig. ^. 

We now turn to studying differential distributions, and start by considering those 
defined in terms of final-state leptons. Observables sensitive to the hadronic activity of 
the events, be either relevant to 6-jets or to i?-hadrons, will follow later. In the case of 
MC-based simulations, several leptons can appear in the final state. We use the MC-truth 
information to select the two which emerge from the hard process, and we shall simply 
refer to them as "the leptons" henceforth. A more realistic analysis may select leptons on 
the basis of their hardness, measured e.g. with their Pt^s; in practice, for the processes we 
are considering here (and thanks to the fact that the 6-hadrons have been set stable) the 
two approaches are equivalent. 

In fig. ^ the invariant mass of the final-state lepton pairs is shown. The effect of the 7* 
contribution to l^i^bb production is clearly visible at small invariant masses. In this plot, 
with limit ourselves to presenting only the aMC@NLO results, since the other simulations 
give results which are essentially identical to the present ones. 

In figs. ^ and ^ the transverse momenta and the pseudorapidities of the charged leptons 
are shown separately according to their electric charges. In the two upper insets we have 
used the same patterns and conventions as in fig. ^ - these will be used throughout this 
paper. In the case of Wbb production, the effects of the NLO corrections are especially 
pronounced at large Pr's, where they are the signal of new partonic subprocesses opening 
up at this order, and in particular of those which include an initial-state gluon, such as 
qg. Results after matching with showers consistently show a similar behaviour. The same 
large enhancement is not present in the case of Z production, which receives gluon-initiated 
contributions already at the LO; again, this trend is seen also after matching with showers. 
The lowest insets (solid magenta curves) show the ratios of the aMC@NLO results relevant 

*We call aMC@LO the analogue of aMC@NLO, in which the short-distance cross sections are computed 
at the LO rather than at the NLO. Its results are therefore equivalent to those one would obtain by using, 
e.g., MadGraph/MadEvent (ETl interfaced to showers. 
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Figure 3: Invariant mass distribution of final-state Icpton pairs, as predicted by aMC@NLO. 
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Figure 4: Transverse momentum of the charged leptons in (I'bb (left panel) and l^l~hh (right 
panel) production, shown separately for positive and negative charges. The upper and middle 
insets follow the same patterns as those in fig. ||. The lower inset (magenta solid histogram) is 
the ratio of the aMC@NLO results relevant to positively-charged leptons over those relevant to 
negatively-charged ones. 



to positively-charged leptons over those relevant to negatively-charged ones. In the case of 
Wbh production the behaviour is similar to what has been described recently in ref. [ p2| , 
while for Zbb this distribution is flat, as expected. 

A complementary aspect of the different parton luminosities that contribute to Wbh 
and Zhb production can be appreciated by looking for example at the transverse momentum 
distributions of the Iv and pairs (i.e. of the virtual W and Z bosons respectively), 

shown in the left panel of fig. ^. In the case of (.ubb final states, the aMC@LO and LO 
results are very close to each other, which is not the case for i'^i'bb production. This is 
due to the fact that the ^(/-initiated channel in the latter case is responsible for much more 
QCD radiation in MC-based simulations than the qq channel (the latter being identical to 
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Figure 5: As in fig. ||, for the pseudorapidity of tlie cfiarged leptons. 

the mechanism that induces Wbb production) . This larger amount of radiation hardens the 
virtual-Z pr spectrum predicted by aMC@LO, making it more similar to the aMC@NLO 
result than in the case of Wbb production. The NLO Pt spectra of the virtual vector bosons 
are closer to the aMC@NLO results, because at that order one does get contributions from 
gluon-initiated channels to both ii'bb and i'^£~bb final states, and thus the relative changes 
obtained when matching with parton showers have a milder impact. This is an example 
of the "stabilizing" pattern that one observes when higher-order perturbative results are 
taken into account. On the other hand, the rapidity distributions of the lepton pairs do 
not change significantly under QCD radiation, as is shown in the right panel of fig. |^. 

We remark that we do not find any significant enhancement in the large-p^ tails of 
vector bosons when going from the NLO to the aMC@NLO predictions, at variance with 



the POWHEG Wbb result of ref. [3£] which necessitates an ad-hoc perturbative tuning for 
this reason. It should further be stressed that the almost perfect coincidence between the 
aMC@LO and LO results for may not occur by simply changing the tuning of the 

shower parameters in Herwig: indeed, we have verified that the spectrum predicted by 
Pythia 6 is slightly different w.r.t. the aMC@LO or LO ones. It is also interesting to 
notice that the features discussed here and that affect the low- and intermediate-pr regions 
of lepton pairs are not visible in the case of the individual lepton pj. spectra, fig. ^ We have 
verified that kinematical correlations are such that, for a fixed small or intermediate value 
of Pt{(-^), one integrates over a range that causes the local differences between 

the aMC@LO and LO results for the latter transverse momentum to be averaged out. 

The left panel of fig. ^ presents the cos 9* distribution, computed by separating the 
positively- and negatively-charged lepton contributions. We remind the reader that such 
an observable is defined as the cosine of the angle between the chosen charged lepton, and 
the direction of flight of the parent vector boson, in the rest frame of the latter. Clearly 
visible are the strong angular correlations and charge asymmetry in the ii/bb case. For 
i^i~bb production such correlations and asymmetries, while present, are much milder, and 
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Figure 6: Transverse momentum (left panel) and rapidity (right panel) of the and t'^ pairs 
(i.e. of the virtual W and Z bosons respectively) in Ivhh and £'^£^bb production. The insets follow 
the same patterns as those in fig. ||. 



likely not observable in a real experiment. For both processes, the aMC@NLO results are 
basically identical to those of aMC@LO, NLO, and LO, and thus we refrain from showing 
the latter here. 

In the right panel of fig. 0, where we consider only leptons with positive electric charge 
to be definite, we plot the ratio of the lepton transverse momentum over the same quantity, 
obtained by imposing a phase-space (i.e., flat) decay of the parent vector boson; hence, 
this ratio is a measure of the impact of spin correlations on the inclusive-lepton pj.. We 
see that differences between correlated and uncorrelated decays can be as large as 20%, 
and vary across the kinematical range considered. This confirms that the inclusion of spin- 
correlation effects is necessary when an accurate description of the production process is 
required. We stress again that our computations feature spin correlations exactly at the 
matrix-element level, including one-loop ones. It is interesting to observe that, while in the 
case of Zbb production all four calculations give similar results (see the lower inset), this 
happens in Wbb production only for Pt{^~^) ^ 50 GeV (see the upper inset). At Pt values 
larger than this, aMC@NLO and NLO predict ratios that differ from the corresponding 
aMC@LO and LO ones. Once again, this is a manifestation of the significant impact of 
gluon-initiated, NLO partonic processes on Wbb cross sections, and is consistent with the 
findings of ref. [^] (relevant to the associated production of W bosons with light jets). 

In figs. P and ^ the transverse momenta and the pseudorapidities of the two hardest 
6-jets are shown. Differences in normalisation are consistent with what we expect on the 
basis of inclusive K factors; differences in shapes are typically small, but visible. We point 
out that for an event to contribute to the hardest-6-jet observables shown here it is sufficient 
that one b-jet be present in the event; the other b quark emerging from the hard process 
can have arbitrarily small momentum. 

In the left panel of fig. [l^, the AR separation between the two hardest 6-hadrons 
(for the MC-based simulations) or between the b and b quarks (for the NLO and LO 
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cosfl- PtD*] [GbV] 

Figure 7: Left panel: cos 0* distribution of final-state charged leptons for different charges. All 
histograms have been obtained with aMC@NLO. See the text for the observable definition. Right 
panel: ratio of the results for the Pt of the positively-charged lepton over the same quantity com- 
puted by neglecting production spin correlations. The insets follow the same patterns as those in 
fig- i 




PT[Jb,i] [GeV] p,[Jt,j] [GeV] 



Figure 8: Transverse momentum of the hardest (left panel) and second-hardest 6-jet (right panel) 
in Wbb and Zbb production. The insets follow the same patterns as those in fig. ^. 

computations) is shown. Differences between the Whh and Zbb processes are manifest. In 
the former case the two Vs originate from a final-state gluon sphtting, and they will thus 
tend to be quite close in pseudorapidity. On the other hand, the two Vs in Zhb production 
can arise from the uncorrelated branchings of the initial-state gluons in the gg channel, and 
in this way they will naturally acquire a large separation in pseudorapidity, which is directly 
related with large-Ai? values. However, a hh pair arising from a final-state gluon branching 
can be easily separated in pseudorapidity by QCD radiation. This is the reason why the 
parton-level LO result in the case of Whh production is so different from the other three 
predictions (as shown by the symbols in the upper inset). Both parton-level NLO (through 
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radiation present at the matrix-element level) and aMC@LO (through radiation due to 
parton showers) results are in fact much closer to aMC@NLO than parton-level LO result 
is. This does not happen in the case of Zbb production, since as discussed before the b and 
b quarks can be well-separated in pseudorapidity already at the LO. It should be stressed 
that the 6-hadrons that contribute to the AR separation shown in fig. |l^ are not subject 
to any lower cuts in p^. Thus, one expects that the effects of extra radiation be diminished 
when imposing a pr cut or, which is equivalent, by studying the same distribution in the 
case of 6-jets. We have verified that this is indeed the case, i.e. that when a minimum-p^ cut 
is imposed on the two 6-hadrons the pattern of NLO QCD corrections in Wbb production is 
more similar to that observed in Zbb production. This is another example of the possibility 
of testing detailed properties of QCD radiation by considering low-p^^ events. It should be 
clear that from the theoretical viewpoint such studies can be sensibly performed only by 
retaining the full 6-mass dependence. 

The right panel of fig. |l^ shows the mass of the 6-jets in the events. The observable 
is inclusive over all 6-jets, which implies that a given event may enter more than once 
in the plot. Striking is the onset of the 66-jet contribution in the Wbb result around 
nij^ w 12 GeV. In the case of Zbb production this effect is almost invisible, a consequence 
of the fact that the fraction of events where a 6-jet is actually a 66-jet is much smaller than 
for Wbb production, see fig. ^. The distribution discussed here measures the activity inside 
a jet, and one cannot expect fixed-order parton-level results, where a jet consists of one 
or two particles, to be particularly sensible in this case. In fact, we see that fixed-order 
results are very different from MC-based ones. On the other hand, the differences in shape 
when going from aMC@LO to aMC@NLO are small, in particular for Zbb production, as 
expected for observables which are insensitive to emissions at large relative Pt's. We also 
point out that the knee at rrij^ ^ 12 GeV would appear as a feature of gluon iets if the 6-jet 
definition of ref. |5C] were used. This stresses again the fact that, at small and moderate 
Pt's, the usual definition gives more intuitive results. On the other hand, at large jet p^'s 
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Figure 10: Left panel: AR separation between the two hardest 6-hadrons (aMC@NLO and 
aMC@LO) or the b and b quarks (NLO and LO) in the event. Right paneh invariant mass of the 
6-jets, inclusive over all 6-jets in the event. The insets follow the same patterns as those in fig. 0. 
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Figure 11: Transverse momentum fraction carried by 6-jets. See the text for details. 



the onset of the 66-jet contribution to nij. is largely smeared out. 



In fig. 11 we show the ratio of the total transverse momentum frbfe] of 6-jets, over the 



total transverse hadronic momentum Pt . In the context of parton-level computations, by 
"hadrons" we simply understand QCD partons. At the parton-level LO, the configurations 
with one 66-jet or with two 6-jets (each of which contains one b quark) give contribution 
at PrUbl/Pr = 1- Configurations with one 6-jet that contains only one b quark contribute 
to 0.5 < P-rijbl/PT < 1 if the other b quark has Pt < 20 GeV (i.e., it is softer than a jet 
is required to be), while values PrUbj/PT < 0.5 can be obtained when the other b quark 
has Pt > 20 GeV and \ri{b)\ > 2.5 (i.e., it is outside the 6-jet tagging region in pseudora- 
pidity). What was said above implies that PrUbl/Pr = 1 is an infrared-sensitive region, 
which gives rise to Sudakov logarithms at higher order; this explains the behaviour of the 



^We stress that Pt is defined without including the underlying event and pile-up contributions. 
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Figure 12: Invariant mass of the pair of the two leading 6-jets. WH{-^ ii^bb), ZH{-^ £^£^bb), 
tvbb, and £^£~bb results are shown, with the former two rescaled by a factor of ten. 

parton-level NLO results there. Furthermore, the LO contributions to the Prijbl/PT < 0.5 
region decrease when increasing the maximum-pseudorapidity cut on jets. This is only 
marginally the case at the NLO (because of the presence of a hard light parton in the real- 
emission contributions), which explains the longer tail of the latter results w.r.t. the LO 
ones. The arguments above obviously do not apply to the context of an event generator; 
this is confirmed by the similarity of the aMC@NLO and aMC@LO results. Firstly, at 
Prijbl/PT = 1 Sudakov logarithms are properly resummed. Secondly, the extra radiation 
generated by parton showers implies that quite a few hadrons will lie outside 6-jets, hence 
shifting further the P-rUbl/PT results to the left of those relevant to parton-level NLO com- 
putations. This shift is also present when passing from the aMC@LO to the aMC@NLO 
predictions in Wbb production, while in the case of Zbb production these two results are 
very similar (up to an overall rescaling by the inclusive K factor) . We are finding here the 
same pattern already discussed for a few observables in this paper. Namely, the opening of 
gluon-initiated partonic channels at the NLO in Wbb production implies a richer hadronic 
activity w.r.t. the corresponding LO case, which is only marginal in the case of Zbb produc- 
tion owing to the dominance of the gg channel already at the LO there. Hence, the relative 
enhancement of the hadronic activity outside the 6-jets when going from aMC@LO to 
aMC@NLO is stronger for Wbb production than is for Zbb production. 

Finally, as a simple application to Higgs searches of the calculations presented in this 
paper, we show in fig. |l^ the invariant mass of the two leading 6-jets in WH{^ ivbb), 
ZH{-^ i'^i'bb), ivbb, and £'^£~bb events. The former two processes (the "signal") have 
been simulated with MC@NLO , with a Higgs mass niH = 120 GeV. The tail at 

®In the process of validating aMCQNLO, we had checked that it gave results identical to MC@NLO for 
all the processes implemented in the latter. Hence, we could have equally well employed aMCQNLO to 
simulate the signal here. 
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m[jh^i, jb,2] > ''T^H is due to the fact that the jet momenta are typically larger than those of 
the 6-hadrons they contain, owing to the contributions of other final-state hadrons emerging 
from initial-state showers. This is compensated by the fact that the 6-hadron momenta 
are only a fraction of those of their parent b quarks, the complementary fraction being 
lost to radiation which may end up outside the jets. These two effects smear the Higgs 
peak. Furthermore, in some events the b quarks entering the two hardest 6-jets do not arise 
from the Higgs decay, but from a g ^ bb branching in the shower phase. Although rare 
indeed, these events may result in invariant masses much larger than the Higgs pole mass. 
The comparison given here is just an example of an analysis in which both the signal and 
its irreducible backgrounds can be computed at the same precision with (a)MC@NLO, 
improving upon both fixed-order and LO-based Monte Carlo descriptions. 



3. Conclusions and outlook 

In this work we have presented results for the ivbb and i^i^bb production processes, 
accurate to the NLO in QCD and that include the matching to parton showers according 
to the MC@NLO formalism. Our approach is fully general, completely automated, and 
opens the way to performing comparisons with experimental data from the Tevatron and 
the LHC at the highest theoretical accuracy attainable nowadays. 

By studying a limited but representative set of observables, we have shown that several 
are the elements to be kept into account in order to achieve reliable and flexible predictions 
for this class of processes: spin correlations of the final state leptons emerging from the 
decays of the vector bosons, heavy-quark mass effects, and a realistic description of the 
final states, obtained thanks to the interface with a shower and hadronisation program. As 
we have seen, NLO QCD corrections have a highly non-trivial impact, since they lead not 
only to large enhancements of total rates, but also to significant changes in the shapes of 
distributions. In this respect, the opening at the NLO of new partonic channels, and in 
particular of those involving gluons, plays a fundamental role. In general and apart from 
well-understood cases in which pure perturbative results are not meaningful, one observes 
that at the NLO level fixed-order and MC-based results are closer to each other than the 
corresponding LO ones. This is in keeping with naive expectations based on perturbation 
theory, and it is significant in that it shows that the very large corrections affecting the 
processes considered here do not pose problems when the matching with parton shower 
Monte Carlos is carried out according to the MC@NLO method. 

Thanks to the aMC@NLO implementation, several QCD issues interesting on their 
own can now be addressed. One example over all is the study of NLO corrections, mass 
effects and radiation pattern in final-state gluon splitting, for which Wbb production offers 
a particularly clean environment. Gluon splitting in the initial state and the role of the b 
PDF (and therefore of different schemes for the predictions of total and differential cross 
sections) can be assessed by considering Zbb production. The outcome of such a study 
can then be applied to the Hbb case. In particular, available predictions in the five-flavour 
scheme at the NNLO for the fully-inclusive production of a Z in association with bottom 
quarks [53|, and for Z+1 6-jet at the NLO |2^, can now compared with our four-flavour- 
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scheme results. In addition, QCD radiation effects on high-p^^ bb pairs, which can be merged 
into one jet, are also of interest in boosted-Higgs searches [54|. Finally, spin correlation 
effects may also be investigated to gather more insight on the production mechanisms 
in QCD, and possibly to distinguish them from other competing hard reactions, such as 
double-parton scatterings. We plan to address some of the above issues in detail in the 
near future. 

We conclude by pointing out that event files relevant to the processes studied in this 
paper (as well as to others) are publicly available at http://amcatnlo.cern.ch. Work to 
make the use of aMC@NLO public from the same site is in progress. 
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